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ABSTRACT 

The objective of this paper is to study an unsteady Magneto hydrodynamic (MHD) viscous, incompressible free 
convective flow of an electrically-conducting, Newtonian fluid through porous medium bounded by an infinite vertical 
porous plate in the presence of radiation, chemical reaction, heat source and surface temperature oscillation and as well as 
concentration oscillation. The governing equations are solved analytically using complex variables. Solutions for the 
velocity field, the temperature field and concentration field are obtained and discussed through graphs. The expressions for 
the Skin friction, the rate of heat transfer in the form of Nusselt number and the rate of mass transfer in the form of 
Sherwood number are also derived and discussed. It is interesting to notice that Steady mean flow velocity is decreased 
with increasing of thermal radiation parameter, magnetic parameter, chemical reaction parameter, Heat source parameter 
and Schmitt number, while it is increased with the increase of Grash of number and modified Grash of number. Unsteady 
main flow velocity is decreased with increasing of thermal radiation, magnetic parameter, chemical reaction, Heat source 
parameter, frequency of oscillation and Schmitt number, while it is increased with the increase of Grash of number and 
modified Grash of number. 
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INTRODUCTION 

Free convection flows are of great interest in a number of industrial applications such as fiber and granular 
insulation, geothermal systems etc. Buoyancy is also of importance in an environment where differences between land and 
air temperatures can give rise to complicated flow patterns. Magnetohydrodynamic has attracted the attention of a large 
number of scholars due to its diverse applications. Convective heat transfer through porous media has been a subject of 
great interest for the last three decades. Kim et al. [1] and Harris et al. [2] solved the problem of natural convection flow 
through porous medium past a plate. Raptis and Kafousias [3] analyzed the effects of magnetic field on steady free 
convection flow through a porous medium bounded by an infinite vertical plate. Raptis [4] investigated time varying 
two-dimensional natural convection flow of a viscous, incompressible and electrically conducting fluid past an infinite 
vertical porous plate embedded in a porous medium. Chamkha [5] studied MHD free convection from a vertical plate 
embedded in a thermally stratified porous medium. Chamkha [6] also investigated MHD free convection from vertical 
porous plate embedded in a porous medium with Hall effects. 
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Verma[7] studied the transient natural convection boundary layer flow, over a semi-infinite horizontal plate with 
the plate temperature oscillating with a non-zero mean, and a stationary free stream using matched asymptotic expansions 
and identifying the presence of four discrete oscillating layers at large distances from the leading edge of the plate. 
Singh et al. [8] have investigated the effects of surface temperature oscillation on transient two-dimensional free 
convective flow through a porous medium bounded by an infinite vertical plate using double asymptotic expansions in 
powers of the frequency parameter. Agrawal et al. [9] studied the influence of free convection and mass transfer on 
hydromagnetic flow past a vibrating infinite isothermal and constant heat flux vertical circular cylinder using Laplace 
transforms for both cooling (Grashof number greater than zero) and heating (Grashof number less than zero) of the 
cylinder). Takhar and Ram [10] considered Hall current and surface temperature oscillation effects on natural convection 
magnetohydrodynamic heat-generating flow. Soundalgekar et al. [11] obtained, exact solutions for heat transfer from an 
infinite vertical oscillating plate with foreign mass injection using the Laplace-transform technique with a linear temporal 
variation of plate temperature showing that Skin-friction increases with increasing frequency, Schmidt number and thermal 
Grashof number, but is reduced with a rise in species Grashof number. 

Das et al. [12] have also investigated unsteady natural convection from a vertical surface with surface temperature 
oscillation effects. Merkin and Pop [13] investigated the unsteady free convection boundary-layer flow near the stagnation 
point of a two-dimensional cylindrical surface embedded in a fluid- saturated porous medium with surface temperature 
oscillation using both numerical and asymptotic solutions showing that external to the thin boundary layer on the surface a 
steady flow is induced for large times and rapid oscillations. Chaudhary and Jain [14] have discussed the combined effects 
of oscillating surface temperature and fluctuating surface velocity, on hydro magnetic convection with internal heat 
absorption. In all the above studies thermal radiation heat transfer was neglected, despite the important contribution this 
mode of heat transfer exerts in high temperature manufacturing processes, fire dynamics phenomena, rocket exhaust 
plumes etc. Anderson et al ., (15) have studied the diffusion of a chemically reactive species from a linearly stretching 
sheet. Anjalidevi and Kandasamy (16) investigated the effect of a chemical reaction on the flow along a semi-infinite 
horizontal plate in the presence of heat transfer. Anjalidevi and Kandasamy (17) have studied the effect of a chemical 
reaction on the flow in the presence of heat transfer and magnetic field. Muthucumaraswamy and Ganesan (18) have 
analyzed the effect of a chemical reaction on the unsteady flow past an impulsively started semi-infinite vertical plate, 
which is subject to uniform heat flux. McLeod and Rajagopal (19) have investigated the uniqueness of the flow of a Navier 
Stokes fluid due to a linear stretching boundary. Raptis et al ., (20) have studied the viscous flow over a non-linearly 
stretching sheet in the presence of a chemical reaction and magnetic field. In 1961, Sakiadis (21) who developed a 
numerical solution for the boundary layer flow field over a continuous solid surface moving with constant speed. 

The aim of the present investigation is to study hydromagnetic free convection of an optically-thin gray gas from 
vertical flat plate subject to a surface temperature oscillation with significant thermal radiation, Chemical reaction and Heat 
source. The fluid is considered to be a gray, absorbing-emitting but non-scattering medium. The effects of the magnetic 
parameter (Hartmann number, M 2 ), free convection parameter (Grashof number, Gr), frequency of oscillation (&>) and the 

radiation parameter (K } ) on primary and secondary velocity distributions, Shear stresses, Nusselt number and Sherwood 
number at the plates are presented graphically and in tabular form. 
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MATHEMATICAL MODEL AND COMPLEX VARIABLE SOLUTIONS 



We consider a two-dimensional free convection effects on unsteady MHD heat and mass transfer characteristics 
of a radiated vertical porous plat moving with constant velocity, U, in the presence of transverse magnetic field. 
The surface temperature of the plate oscillates with small amplitude about a non-uniform mean temperature and the 
concentration oscillation is also considered. The co-ordinate system is such that the x*- axis is taken along the plate and 
y* - axis is normal to the plate. All the fluid properties are considered to be constant except the influence of the density 
variation in the buoyancy term, according to the classical Boussinesq approximation. Since magnetic Reynolds number is 
very small for metallic liquid or partially ionized fluid the induced magnetic field produced by the electrically-conducting 
fluid is negligible. The Hall effects, the viscous dissipation and the joule heating terms are also neglected. The governing 
equation for this investigation is based on the balances of mass, linear momentum, energy, and species concentration. The 
equations can be written in Cartesian frame of reference, as follows. 
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Where w* and w* are x*-direction and z*-direction velocity components, p is the fluid density, p is the viscosity, c p 
is the specific heat at constant pressure, a is the fluid electrical conductivity, P is the volumetric coefficient of the thermal 
expansion, pi is the volumetric coefficient of the concentration expansion, B 0 is the magnetic induction, K* is the 
permeability of the porous medium, T * is the dimensional temperature, D 7 is the coefficient of chemical molecular 
diffusivity, C* is the dimensional concentration, k is the thermal conductivity of the fluid, g is the acceleration due to 
gravity, ^ is the local radiative heat flux, k £ the reaction rate constant respectively. 



The term Q 0 (T -TJ is assumed to be amount of heat generated or absorbed per unit volume go is a constant, 
which may take on either positive or negative values. When the wall temperature T * exceeds the free stream 
temperature 1^, the source term Q 0 > 0 and heat sink when Q 0 < 0. 



The boundary conditions for the velocity, temperature, and concentration fields are given as follows 
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Where (f) is designated as wall-free stream temperature difference, %j) is designated as free stream concentration 

U 

difference, c =— dimensionless velocity ratio and 0} is the frequency of oscillation in the surface temperature of the 

S 



plate. 



The conditions (6) suggest solutions to equations (1) to (5) for the variables U* p W * f (j) and 'ifj of the form. 

it' = + ee iw ’*'i4 

w* = wj + ee lw '*V; 

<j, = T‘-TZ = 
t = c* - q = e„<;c$ + 

For the case of an optically-thin gray gas, the thermal radiation flux gradient may be 
Expressed as follows (Siegel and Howell 1993): 



(7) 



3y : 



= 4 au* (C - r 4 ) 



( 8 ) 



And rj * is the radiative heat flux, a is the absorption coefficient of the fluid and C7* is the Stefan-Boltzmann 
constant. We assume that the temperature differences within the flow are sufficiently small such that 7 " * 4 may be 

expressed as a linear function of the temperature. This is accomplished by expanding T*^m a Taylor series about T^ 4 
and neglecting higher order terms, leading to: 

r 1 = 4 r; 3 r - sr; 4 « 

|2( = 16affT«<T*-ri) do) 

The components Hq, represent the steady mean flow, represents temperature fields, represents 
concentration fields and satisfy the following equations: 
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The corresponding boundary conditions are 

K=u, w* = o, t* = t:,c* = Q , 
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The components Wj represent the unsteady mean flow, £7^ represents temperature fields, represents 



concentration fields and satisfy the following equations: 
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The corresponding boundary conditions are 

ui = u, Wj = o, r = tz , c 
t£->o, Wi -> o, r -ir^c ->c ( 

We introduce following dimensionless quantities to normalize the flow mode 
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Where G r is the Grash of number, C m is modified Grash of number, M is Magnetic parameter, K x is the thermal 

radiation-conduction number, is the chemical reaction parameter, S-l is dimensionless temperature and Ci is 
dimensionless concentration. Using equation (23), the dimensionless form of the equations are given by 
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The corresponding boundary conditions for steady mean flow are 

u Q = t, w Q = 0, & Q = 1, C 0 = 1 at y = 0 

% = Qj Wq = 0, 0q, = 0, C Q = 0, at y -? oo (32) 

The corresponding boundary conditions for unsteady mean flow are 
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it 1 = 0,w 1 = Q, 9 1 = 0, C 1 = 0, at y -*■ ot> 

Introducing the complex variable = + uVq, ti = + iw 1 , eqs (24), (25), (28) and (29) implies 

f= -6 r e,-c„c 0 

0-(to + M 1 )H-C r 8 l -6 r C 1 

The corresponding boundary conditions are: 

F = 1, Sq = 1, Cq = 1 at y = 0 

f = 0 , 8 q = 0, C(j = 0, at y -» oo 
J/ = 1, Sj = 1, Cj = 1 at y = 0 

H = 0 , Sj = 0, C 1 = 0, at y -> oo 
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(33) 



(34) 



(35) 
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Solutions 

Solving the above set of equations with respect to the corresponding boundary conditions, the following solutions 
are obtained. 
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Equation (34) subject to the boundary conditions (36) can be solved and solution for the steady mean flow is be 
expressed as 
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Wi(y,t) = [m 9 SinB^' — m^osB^yle A±y + m 7 SinB 2 y + m s CosB 2 y]e A * y + 

[m 5 SinB 3 y+ m 4 CosB^y]s~ Aiy 

The functions Sq. and 8 ^ denote the temperature fields due to the main flow and cross flows, respectively. 



(45) 



Of interest in practical MHD plasma energy generator design are the dimensionless shear stresses at the plate, 
which may be defined for steady and unsteady mean flow, respectively, as follows: 
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It is evident from equations (46) and (47) that the shear stress component due to the main flow for the steady 
mean flow given by equation (46) and the shear stress components due to main and cross flows given by equation (47) do 
not vanish at the plate. Inspection of these expressions also reveals that the shear stress component as defined by equation 
(46) due to steady mean flow is subjected to a non-periodic oscillation that depends on Hartmann number and radiation 
conduction parameter. In contrast to this, the shear stress components as computed in equation (47) due to the main and 
cross flows for an unsteady mean flow are subjected to periodic oscillation which is a function of not only Hartmann 
number and radiation-conduction parameter, but also the Prandtl number and the frequency of oscillation. The shear stress 
for the equation (46) will vanish at the plate (y=0) at a critical value of the free convection parameter i.e. Grash of number, 
defined by the conduction: 



Cr crit = v M i(v 5 i + Wi) And ° 

The shear stress for equation (47) will vanish at the plate (y=0) When: 
= ( A i + +- A i ) + *((^1 + B 2 )] and G m = 0 



(48) 



(49) 



Also of interest in plasma MHD generator design is the dimensionless temperature gradient at the plate. This can 
be shown to take the form, for the unsteady main flow, as follows: 






(^)y=0 = "A = V* + *1 



For the unsteady cross flow the dimensionless temperature gradient at the plate (y=0) is: 



(50) 
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For the unsteady main flow and cross flow the dimensionless concentration gradient at the plate (y=0) is given by: 
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(52) 
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(~£)y=o = -V 5 4 = - V 5 2 + ^ (53) 

RESULTS AND DISCUSSIONS 

The selected computations for the velocity, temperature and concentration for different values of the material 
parameters have been presented in figures 1 - 20. Computations for the shear stresses, temperature gradient at the plate for 
steady and unsteady mean flow are provided in Tables 1-5. We note that in figures 1-6, only steady mean flow is 
simulated, for which there will be no surface temperature oscillation. Figures 7-14, correspond to the unsteady mean flow 
distribution due to the main flow (ui) and figures 15-20 corresponds to the unsteady mean flow due to cross flow (wi). 

An increase in radiation conduction number Ki has an adverse effect on the velocity due to steady mean flow (u 0 ) 
for all values of y, as shown in figure 1 . represents the relative contribution of thermal radiation heat transfer to thermal 
conduction heat transfer. For K x < 1 thermal conduction exceeds thermal radiation and for K x > 1 this situation is reversed. 
For Ki = 1 the contribution from both modes is equal. From figure 2, it is observed that, an increase in free convection 
parameters, Grash of number (Gr), increases the velocity due to steady mean flow. Similar effect has been observed from 
figure 3, in the case of modified Grash of number Gm. An increase in the magnetic parameter (M) has an inhibiting effect 
on steady mean flow as shown in figure 4. From figure 5, it is noticed that, as chemical reaction parameter ( )') 

increases, the velocity due to steady mean flow tends to decrease. An increase in Schmitt number (sc), has an adverse 
effect on the velocity due to steady mean flow (u 0 ) as shown in figure 6. 

From figure 7, it is observed that, an increase in radiation conduction number Ki decreases the unsteady mean 
flow (ui). Whereas cross flow velocity (w0 increases as increases (as shown in Figure 15). The effects of free 
convection parameters, Grash of number (Gr) and Modified Grash off number on unsteady mean flow are similar to for 
steady mean flow as shown in figures 8 and 9. An increase in the magnetic parameter M, has an adverse effect on the 
velocity due to unsteady mean flow as shown in figure 10. From figure 11, it is noticed that, as chemical reaction 
parameter ( y) increases, the velocity due to unsteady mean flow tends to decrease. From figure 12, it is observed that, 

as Prandtl number ( Pr) increases, the velocity due to unsteady mean flow tends to decrease. The effect of oscillation 
frequency (oj) is computed, in figure 13, on the unsteady mean flow velocity component decreases as increases. From 
figure 14, it is observed that, an increase in Schmitt number (sc), decreases the velocity due to unsteady mean flow. 

From figure 15, it is noticed that, an increase in Grash off number, decreases the velocity due to unsteady mean 
flow distribution due to cross flow (wi), where as reverse effect is observed from figure 16 as modified Grash off number 
increases. From figure 17, unsteady mean cross flow velocity decreases in magnitude with a rise in magnetic parameter M. 
Unsteady mean flow distribution due to cross flow decreases as oscillation frequency (*J) increases as shown in figure 18. 

The effects of M, K b Gr, Gm, Sc, &> on share stress at the plate due to steady mean flow, unsteady main flow and unsteady 

cross flow are numerically shown in Table 1, Table 2 and Table 3 respectively. Also the effects of K b and S on the 
temperature gradient at the plate due to unsteady main flow are numerically shown in Table 4. 

From Table 1 it is noticed that an increase in M, Ki and Gm results in a decrease in share stress at the plate due to 
steady mean flow, while it increases as an increase in Gr and Sc. From Table 2, it is observed that an increase in M, KI, 
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Gr, Sc and OJ leads to an increase in magnitude in the share stress at the plate due to unsteady main flow, while reverse 

effect is noticed for an increase in Gm. From Table 3, it is noticed that an increase in M and Ki leads to an increase in 
magnitude in the share stress at the plate due to unsteady cross flow while reverse effect is noticed for an increase in Gr, 
Gm, Sc and (x). From Table 4, it is noticed that an increase in K x and S leads to an increase in the temperature gradient at 
the plate due to unsteady main flow. 





Figure 1: Velocity Distribution Due to a Steady Mean Flow Figure 2: Velocity Distribution due to a Steady 
(u 0 ) for Various Thermal Radiation- Conduction Numbers Mean Flow (uO) for Various Grashof Numbers (Gr) 

(Ki) 




Figure 3: Velocity Distribution Due to a Steady Mean Flow 
(UO) for Various Modified Grash of Numbers (Gm) 




Figure 4: Velocity Distribution Due to a Steady 
Mean Flow (uO) for Various Values of Magnetic 
Parameter (M) 
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Figure 5: Velocity Distribution Due to a Steady 
Mean Flow (uO) for Various Values of Chemical Reaction 

(y) 



Figure 6: Velocity Distribution Due to a Steady 
Mean Flow (uO) For Various Values of Schmitt 
Number (Sc) 




o ■ ~~~ T ~ — 

1 2 : 3 4 5 




Figure 7: Unsteady Mean Flow distribution Due to Main 
Flow (ul) For Various Thermal Radiation- Conduction 
numbers (kl) 



Figure 8: Unsteady Mean Flow Distribution Due to 
Main Flow (ul) For Various Values of Grashof 
numbers (Gr) 
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Figure 9: Unsteady Mean Flow Distribution Due to Main 
Flow (ul) For Various Values of Modified Grashof 
Numbers (Gm). 



Figure 10: Unsteady Mean Flow Distribution Due 
to Main Flow (ul)for Various Values of Magnetic 
Parameter (M) 





Figure 11: Unsteady Mean Flow Distribution Due to Main Figure 12: Unsteady Mean Flow Distribution Due 

Flow (ul) For Various Values of Chemical Reaction (y). to Main Flow (ul) For Various Values of 

Frequency of Oscillation (co). 
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Figure 13: Unsteady Mean Flow Distribution Due to Main Figure 14: Unsteady Mean Flow Distribution Due 
Flow (ul) For Various Values of Prandlt Numbers (Pr) to Main Flow (ul) For Various Values of Schmitt 

Number (Sc) 





Figure 15: Unsteady Mean Flow Distribution due to Cross 
Flow (wl) For Various Thermal Radiation- Conduction 
Numbers (kl) 



Figure 16: Unsteady Mean Flow Distribution Due 
to Cross Flow (wl) For various Values of Grash of 
Numbers (Gr) 
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-3 

X 10 





Figure 17: Unsteady Mean Flow Distribution due to Cross 
Flow (wl) For Various Values Modified of Grash of 
Numbers (Gm) 



Figure 18: Unsteady Mean Flow Distribution Due 
to Cross Flow (wl) For Various Values of Magnetic 
Parameter (M) 




Figure 19: Unsteady Mean Flow Distribution Due to Cross 
Flow (wl) For Various Values of Chemical Reaction (y) 



Figure 20: Unsteady Mean Flow Distribution due to 
Cross Flow (wl) For Various Values of Frequency 
of Oscillation (o) 
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Figure 21: Unsteady Mean Flow Distribution Due to Cross Flow (wl) For various values of S 



Table 1: Shear Stress at the Plate Due to Steady Mean Flow (u 0 ) 
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K, 


Gr 


Gm 


Sc 


f du 0 ^ 
l dy Jy=0 


5.00 


0.50 


7.00 


4.00 


0.50 


1.8730 


7.00 


0.50 


7.00 


4.00 


0.50 


0.9147 


5.00 


1.00 


7.00 


4.00 


0.50 


1.7886 


5.00 


0.50 


9.00 


4.00 


0.50 


3.5656 


5.00 


0.50 


7.00 


5.00 


0.50 


1.5216 


5.00 


0.50 


7.00 


4.00 


1.00 


3.9508 



Table 2: Shear Stress at the Plate Due to Unsteady Mean Flow Due to Main Flow (ui) 



M 


Kj 


Gr 


Gm 


Sc 


CD 


f du ^ 

l dy Jy=0 


5.00 


0.50 


7.00 


4.00 


0.50 


0.50 


-0.1807 


7.00 


0.50 


7.00 


4.00 


0.50 


0.50 


-0.7068 


5.00 


1.00 


7.00 


4.00 


0.50 


0.50 


-0.2019 


5.00 


0.50 


9.00 


4.00 


0.50 


0.50 


0.2545 


5.00 


0.50 


7.00 


5.00 


0.50 


0.50 


0.0553 


5.00 


0.50 


7.00 


4.00 


1.00 


0.50 


-0.3072 


5.00 


0.50 


7.00 


4.00 


0.50 


2.00 


-0.7630 



Table 3: Shear Stress at the Plate Due to Unsteady Mean Flow Due To Cross Flow (wi) 



M 


Kx 


Gr 


Gm 


Sc 


n 


f dwA 

l d y J y=0 


5.00 


0.50 


7.00 


4.00 


0.50 


0.50 


-0.6958 


7.00 


0.50 


7.00 


4.00 


0.50 


0.50 


-0.5313 
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Table 3: ContcL 



5.00 


1.00 


7.00 


4.00 


0.50 


0.50 


- 0.6737 


5.00 


0.50 


9.00 


4.00 


0.50 


0.50 


- 0.8677 


5.00 


0.50 


7.00 


5.00 


0.50 


0.50 


- 0.7430 


5.00 


0.50 


7.00 


4.00 


1.00 


0.50 


- 0.7011 


5.00 


0.50 


7.00 


4.00 


0.50 


2.00 


- 1.0356 



Table 4: Temperature Gradient at the Plate Due to Steady Mean Flow (0 O ) 



Kl 


s 


UeA 

l dy Jy= 0 


0.50 


1.00 


1.2247 


1.00 


1.00 


1.4142 


0.50 


2.00 


1.7321 



CONCLUSIONS 

Exact solutions have been derived using complex variables for the transient magneto hydrodynamic (MHD) 
convection flow of an electrically-conducting, Newtonian, optically-thin fluid from a flat plate with thermal radiation, 
chemical reaction, heat source and surface temperature oscillation effects. Our analysis has shown that: 

• Steady mean flow velocity, u 0 , is decreased with increasing of thermal radiation magnetic parameter (M), 
chemical reaction (y),Heat source parameter (S) and Schmitt number (Sc) , while it is increased with the increase 
of Grashof number (Gr) and modified Grashof number ( Gm ). 

• Unsteady main flow velocity (u^ is decreased with increasing of thermal radiation (K {), magnetic parameter (M), 
chemical reaction (y),Heat source parameter (S), frequency of oscillation (co) and Schmitt number (Sc) , while it is 
increased with the increase of Grashof number (Gr) and modified Grashof number (Gm). 

• Conversely cross flow velocity (wj) is increased with increasing of thermal radiation (Ki ), magnetic parameter 
(M), Heat source parameter (S) and modified Grashof number (Gm), while it is decreased with the increase of 
Grashof number (Gr), chemical reaction parameter (y) and frequency of oscillation (co). 

• An increase in M, radiation-conduction parameter Kj and modified Grashof number (Gm) causes a decrease in 
the shear stress at the plate due to a steady mean flow (u 0 ). 

• An increase in frequency of oscillation (co), magnetic parameter (M), radiation-conduction parameter Kj , Grashof 
number (Gr), Schmitt number (Sc) and modified Grashof number (Gm) increases shear stress at the plate due to 
unsteady mean flow of the main flow ( uj ). 

• Shear stress at the plate with unsteady mean flow due to cross flow (wj) also decreases with an increase of 
magnetic parameter (M ) and radiation-conduction parameter Kj , but it is considerably increased with increase of 
Grashof number (Gr), Schmitt number (Sc), modified Grashof number (Gm) and frequency of oscillation (co). 

• An increase in thermal radiation-conduction parameter (Kj) and Heat Source parameter (S) increases the 
temperature gradient at the plate due to unsteady mean flow (0 O ) . 



Impact Factor (JCC): 4.6257 



NAAS Rating: 3.80 






Analytical Study of MHD Free Convective Heat and Mass Transfer 

Flow Bounded by an Infinite Vertical Plate with Thermal Radiation and Chemical Reaction 



45 



• In the absence of chemical reaction, concentration and heat source parameters, these results are in good agreement 
with the results of O. Anwar Beg and S. K. Ghosh [26]. 

APPENDICES 

= M + ^ , S 1 = S + k L , S 2 = fft , S 3 = S 1 + E£uP f , S 4 = S 2 + icoS e 



a l = (ft - ft = "(ft - 1) , = (ft - *ft). ft = *>(ft - 1), 

A = [(«i ! + + i 

= £ [(V + +5 1 ] Vj . B, = i [C^ 2 + a. z P r 2 ) V2 - SjVa 



^ = J= [(V + + Jj] 1 /! , B, = i [(S 2 2 + ^ 2 - s,] 1 /! , m, - 

lrL2 = ^' m2 = cl+mi+jn2) ' = = 
= = [C S +C 4 ], m s =[l+C 7 + C E ] 



7 
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